Time reversal invariant realization of the Weyl semimetal phase 



Gabor B. Halasz^^^ and Leon Balents^ 
^ Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 
^Trinity College, University of Cambridge, Trinity Street, Cambridge CB2 ITQ, UK 

We propose a realization of the Weyl semimetal pliase tliat is invariant under time reversal and occurs due 
to broken inversion symmetry. We consider both a simple superlattice model and a more realistic tight-binding 
model describing an experimentally reasonable HgTe/CdTe multilayer structure. The two models have the same 
underlying symmetry, therefore their low-energy features are equivalent. We find a Weyl semimetal phase 
between the normal insulator and the topological insulator phases that exists for a finite range of the system 
parameters and exhibits a finite number of Weyl points with robust band touching at the Fermi level. This phase 
is experimentally characterized by a strong conductivity anisotropy and topological surface states. The principal 
conductivities change in a complementary fashion as the system parameters are varied, and the surface states 
only exist in a region of momentum space that is determined by the positions of the Weyl points. 



I. INTRODUCTION 

In the last decades, topological phases of matter have been 
in the focus of intense theoretical and experimental study: for 
a review, see Ref. Q] and references therein. The order ex- 
hibited by these phases is not associated with spontaneous 
symmetry breaking, and it can be described by topological 
invariants that are insensitive to smooth changes in the sys- 
tem parameters;^ As a generic feature, these phases also have 
topologically protected edge states. 

The field of topological phases was revolutionized by the 
discovery of two-dimensional (2D) topological insulators,i^ 
and the subsequent generalization to three-dimensional (3D) 
topological insulators;^ These materials exhibit a bulk energy 
gap between the valence and the conduction bands, similarly 
to normal insulators. On the other hand, they have gapless 
surface states that are topologically protected, therefore con- 
duction is possible on the surface. Since topological insula- 
tors arise due to strong spin-orbit coupling, their prevalence is 
larger within materials consisting of heavier elements;^ They 
find potential applications in the areas of spintronics and quan- 
tum computation. 

It is a recent development that topologically protected sur- 
face states can also be achieved in materials without a bulk en- 
ergy gap: these are the Weyl semimetals.^"- They have band 
touching between the conduction and the valence bands at 
the Fermi level. The band touching points are called Weyl 
points because the dispersion relation around them is linear 
and hence the excitations are equivalent to Weyl fermions. 
Weyl points can have positive or negative helicities, and they 
always appear in pairs. To achieve robust band touching that 
can not be removed by an infinitesimal perturbation, Weyl 
points of opposite helicities must be separated in momentum 
spaccp This requires breaking either the time reversal or the 
inversion symmetry of the system 

Recent papers on Weyl semimetals have predominantly 
studied the case with broken time reversal symmetry;^ One 
notable exception is Ref. HI where time reversal symmetry 
remains intact and inversion symmetry is broken. It was ar- 
gued that a gapless phase appears in three dimensions between 
the normal insulator (NI) and the topological insulator (TI) 
phases. In this paper, we are also interested in the time re- 



versal invariant case, but address specifically how this phase, 
which is in fact the Weyl semimetal, may be designed in a 
NI/TI superlattice. We propose two models: a simple super- 
lattice model adapted from Ref. [H and a more realistic tight- 
binding model describing a HgTe/CdTe multilayer structure. 
The former model is presented in Section and the latter one 
is presented in Section |III] The most prominent experimental 
features are discussed in Section lTVl while the overall conclu- 
sions of the paper are summarized in SectionlV] 



II. SUPERLATTICE MODEL 
A. General description 

The model considered in this section is based on the multi- 
layer structure in Ref. 8: a periodic superlattice of NI and TI 
layers grown in the z direction. It is a simplified tight-binding 
model where we only take the surface states located at the 
NI/TI interfaces into account. These states are labeled by the 
unit cell index and the parallel 2D momentum k = {kx,ky). 

It is known that a realization of the Weyl semimetal phase 
requires breaking either the time reversal or the inversion sym- 
metry of the system. Since we intend to keep the time reversal 
symmetry intact, the inversion symmetry must be broken. To 
achieve that, we introduce a finite voltage V between the top 
and the bottom NI/TI interfaces in each unit cell. The Hamil- 
tonian of the multilayer structure is then 



(1) 



where the Pauli matrices a = (cr^, cr^, cr^) act on the real spin 
degree of freedom and the Pauli matrices f ~ {t^ ^t^ ,t^) 
act on the top/bottom surface pseudospin degree of free- 
dom. The first term describes the NI/TI surface states with 
isotropic Fermi velocity vp, the second term represents the 
inversion-breaking voltage, and the remaining terms describe 
hopping between neighboring interfaces. The hopping ampli- 
tude through a TI layer is At and that through a NI layer is 
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Ajv- In general, both At and An can be functions of the 
parallel momentum k, and the symmetries of these functions 
determine the symmetry of the system. 

The Hamiltonian in Eq. ([T]) can be solved by exploiting the 
translational symmetry in the z direction, and introducing the 
corresponding 3D momentum k = {kx, ky,kz)- By doing so, 
we find that the band dispersion relation is 



El{k) = A^{k,)+[V±VF\k\y 



(2) 



where A{kz) = a/ + + 2 At An cos{kzd) and d is 
the periodicity of the superlattice. The four bands are non- 
degenerate when k 7^ 0, and band touching between the two 
middle bands takes place when = 0. If we assume with- 
out loss of generality that A^ and A jv are both positive, this 
happens when kzd = tt. At = An, and V = u_F|k|. 



B. The Weyl semimetal phase 

If At and Ajv are independent of k, the band touching 
occurs along a circle of radius V/v p in the kz = n/d plane. It 
marks the transition between the NI and the TI phases of the 
material at At = An- We can argue on physical grounds that 
At > Ajv (thin TI layers and thick NI layers) corresponds to 
the NI phase, while At < Ajv (thick TI layers and thin NI 
layers) corresponds to the TI phase. 

However, this band touching is not robust because it re- 
quires the fine-tuning of the condition At = An- To achieve 
robust band touching, we need to make the hopping ampli- 
tudes depend on the momentum k: 
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Furthermore, we can not keep the continuous rotational sym- 
metry around the z axis because then At and An are still 
constants at |k| = V/vf, the only region where band touch- 
ing is possible. On the other hand, the continuous rotational 
symmetry is broken in real crystals as well, and one is only 
left with a discrete rotational symmetry. In the following, we 
demonstrate robust band touching in the reasonable cases of 
the four-fold and two-fold rotational symmetries. 



C. Four-fold rotational symmetry 

In the first case, we assume a four-fold rotational symmetry 
around the z axis and four planes of reflection symmetry: the 
{.T, z} plane, the {j/, z} plane, and the two planes halfway in 
between. These are the natural symmetries of many real ma- 
terials with tetragonal crystal structures. By neglecting any 
contributions depending on |k| only, the lowest order term 
having all the above symmetries and time reversal symme- 
try is cx {k^ + ky)- The k dependent parts of the hopping 
amplitudes are then 



where the polar coordinates k^ = |k| cos and ky = |k| sin 9 
are introduced. The difference At — An depends on the angle 
6 at |k| = V/vf, therefore band touching with At = Ajv 
only occurs at specific points of the circle. The band touching 

also becomes robust because the parameters A^^^ and St,n 
do not require fine-tuning: an infinitesimal change in any of 
them only gives an infinitesimal change in 6, displacing the 
band touching points along the circle. 

Contrary to the case with continuous rotational symmetry, 
now there is a Weyl semimetal phase between the NI and the 
TI phases that exists for a finite range of the parameter val- 
ues. This phase features a finite number of Weyl points at 
which band touching between the two middle bands occurs. 
To be more precise, the solution of the band touching equa- 
tion At = Ajv for the angle 9 is 
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This expression gives 8 Weyl points which are related to each 
other by the symmetry transformations of the system. Since 
I cos(4 6)\ < 1, the condition for the Weyl semimetal phase 

becomes 1/2 < (A^^ - A^°^) / [|k|'^((5T - Sn)] < 1, where 
|k| = V/vf as always in this subsection. Let us assume with- 
out loss of generality that St > Sn, and imagine decreas- 
ing A^^ gradually while keeping the other parameters con- 
stant. This corresponds to a transition from the NI phase to 
the TI phase. The Weyl points then first appear at the lines 
kx = zkky, move along the circle of radius V/vf, and finally 
disappear at the lines k^ = and ky = 0. For an illustration 
of this, see the top half of Fig. [T] 

It can be verified that the band touching points occurring 
in this scenario are indeed proper Weyl points around which 
the band dispersion is linear in all directions. To obtain physi- 
cally transparent results, we assume that St—Sn is sufficiently 
small so that At — A^v is almost independent of 9- This differ- 
ence is then only important along circles of constant kz and |k| 
where there would be no difference otherwise. Consequently, 
the principal directions are the axial (z), the radial (?■), and 
the tangential (t) directions, and an expansion of around a 
band touching point reads 
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A^^]v(k) = StMM^ (cos'' 9 + sin^ 



(4) 



where the effective Fermi velocities corresponding to the prin- 
cipal directions are Vz = d^/ AtAn = ^At, Vr = vf, and 
Vt ~ |kp((5T — Sn) \ sin(4 0)|. The expression in Eq. ^ in- 
deed gives a linear band dispersion in all directions. We can 
now establish that the approximation of small St — Sn re- 
quires Vt <C Vr, i.e. V^{St — Sn) ^ Vp. This is satisfied 
in the reasonable case when the inversion-breaking voltage V 
and the coefficients St,n are small. 

We note that the Weyl points related to each other by rota- 
tions have identical helicities, while those related to each other 
by reflections have opposite helicities: this implies that there 
are 4 Weyl points of each helicity. If we pair up all the Weyl 
points into pairs of opposite helicities, the sum of the result- 
ing separation vectors is zero. This property follows from the 
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FIG. 1 : (Color online) Illustration of the Weyl semimetal phase in the 
cases of four-fold (top) and two-fold (bottom) rotational symmetries. 
We set kz = vr/d in all subfigures. (a) Arrangement of the Weyl 
points with positive (red plus) and negative (blue minus) helicities. 
(b) Trajectories of the Weyl points as the transition from the NI phase 
to the TI phase takes place. The Weyl points first appear at the red 
squares and finally disappear at the yellow diamonds. 



general notion of time reversal symmetry which also implies 
that the anomalous Hall conductivity vanishes. 



D. Two-fold rotational symmetry 

In the second case, we have a two-fold rotational symmetry 
around the z axis and two planes of reflection symmetry: the 
{.T, z} plane and the {y, z} plane. This case is particularly 
important for us because the tight-binding model described in 
Section ini] has the same symmetries. The lowest order terms 
obeying all these symmetries and time reversal symmetry are 
(X {k'^ + ky) and cx (fc^ — ky). However, the former one only 
depends on |k|, and hence it would not break the continuous 
rotational symmetry on its own. For the sake of simplicity, we 
consider the special case of 



A^r^^^ik) = 6t^n kl = ^T.ivlkp cos^ 9, 
and obtain that the solution of At = A^r is 
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Now there are 4 Weyl points in the Weyl semimetal phase that 



occurs when < (A)^' - A^")/ [|k|2(5T - 6j^)] < 1. If 
there is a transition from the NI phase to the TI phase due to 
a gradual decrease in A^', the Weyl points first appear at the 
kx = Q line, move along the circle of radius V/vp, and finally 



disappear at the fc^ = line. For an illustration of this, see the 
bottom half of Fig. [T] 

If i5t — (^AT is sufficiently small, the principal directions 
around the Weyl points are the axial, the radial, and the tan- 
gential directions. Eq. (|6]l is therefore valid in this case as 
well, and the effective Fermi velocities in the principal direc- 
tions are given by Vz = d^/ AtAn = dAr, Vr = vp, and 
Vt = |k|(<5T — SN)\sm{26)\. The approximation of small 
St — Sn holds when vt ^ Vr, i.e. when V{6t — Sn) ^ Vp- 



III. REALISTIC TIGHT-BINDING MODEL 
A. Formulation of the model 

In this section, we consider a periodic multilayer structure 
of strained HgTe and CdTe layers which are grown on top of 
each other in the z direction. This model is in fact a concrete 
realization of the superlattice structure described in Section HII 
because CdTe is a NI and HgTe becomes a TI under straiujii 
The band structures of these materials are well known, and can 
be reproduced with high accuracy from realistic tight-binding 
models. Here we adapt the ten-band tight-binding model de- 
scribed in Ref. [l2l which assumes two s orbitals (s, s*) and 
three p orbitals (p^, Py,Pz) on each atom. 

Both HgTe and CdTe have zinc -blende structures: the an- 
ions (Te) form a face-centered cubic lattice, and the cations 
(Hg, Cd) are located at the positions 1, 1]. This implies 
that each anion (cation) is tetrahedrally coordinated by 4 near- 
est neighbor cations (anions). We assume that the thicknesses 
of the HgTe and CdTe layers are A^i and N2 as measured in 
units of the cubic lattice parameter a. When cutting through 
the structure along the z direction, one finds subsequent lay- 
ers of one atomic thickness consisting of only anions and only 
cations, respectively. The anionic layers all contain Te, while 
there are 2Ni cationic layers containing Hg and 2N2 cationic 
layers containing Cd in each superlattice period. 

The atomic orbitals are labeled according to \R,u,t,a), 
where R is the position of the atomic site, u = {Te, Hg, Cd} 
is the type of the atom, t = {s, s* ,px,Py,Pz} is the type of 
the orbital, and cr = {t, i} is the spin quantum number. The 
Hamiltonian of the system can be written as 



H — Ho + Hj + Hso, 



(9) 



where the first term contains the bare energies of the atomic 
orbitals, the second term describes the interaction (hopping) 
between them, and the third term represents spin-orbit cou- 
pling. These terms are 

Ho = J2 \R,<R),t,a)E^^j^^/R,uiR),t,al (10) 

R,t,(7 

R,t,(7 R',t' 

F,o = \R,u{R),t,a)2X^^j^^L ■ c7{R,u{R),t' 

R,t.t'a.a' 
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where the second sum in Hj goes over all 4 nearest neighbors 
R' of the atomic site R, and u{R) denotes the type of the atom 
at the position R. 

The hopping amplitudes Vu,u' ,t.t' between different types 
of p orbitals are related to each other by the geometry of the 
crystal structure. In particular, they are affected by the uni- 
axial strain e which is defined as the relative elongation of 
the lattice constant in the z direction with respect to those in 
the X and y directions. This strain occurs because the sub- 
sequent layers of HgTe and CdTe are grown on top of each 
other, and there is a slight lattice constant mismatch^ Since 
the lattice constant of CdTe is approximately 0.3% larger than 
that of HgTe, and both materials have a Poisson's ratio 0.5, 
we assume that the relationship between the strains in them is 
e(CdTe) = e(HgTe) + 0.009. Even if a small strain does not 
change the distance between neighboring atoms, the direction 
vector connecting them changes, leading to a different overlap 
between any two orbitals if at least one of them is a p orbital. 
Using simple geometry, the different hopping amplitudes in- 
volving p orbitals are then expressed as 



to each other at this point. The spin-orbit coupling strength 
Atc is averaged in the same way, and the final values of the 
tight-binding parameters are presented in TablesHlandHIl 
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Note that the subscripts u and u' are suppressed for the sake 
of compactness, the label s can stand for both s and s* , and 
all terms are expanded up to first order in e. 
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Te 


Hg 


Cd 




-9.75 


-1.40 


-1.42 




0.12 


4.30 


3.48 




6.08 


6.50 


6.67 




0.333 


0.286 


0.013 



TABLE I: Bare energies Eu,t and spin-orbit coupling strengths An 
for different atom types (all numbers are in eV units). 

The concrete tight-binding parameters are based on those 
in Ref. [13, but they are normalized according to a consistent 
procedure. The bare energies of all orbitals in the CdTe model 
are first shifted such that the Te orbitals have the same aver- 
age energy in HgTe and CdTe. This corresponds to matching 
the arbitrary zero energy levels of the independent HgTe and 
CdTe tight-binding models. The bare energies of the respec- 
tive Te orbitals are then obtained by averaging those in the 
normalized HgTe and CdTe models which are already close 



u/u' 



Te/Hg 



Te/Cd 



u,u' ,p,s,a 



V, 



K 



u.u ,p,s ,a 



u,u' jp,p,a 
^^Uju' ,p.p,7r 



-0.817 
1.044 

-1.404 
1.524 

-0.140 
2.180 

-0.549 



-1.195 
0.753 

-2.064 
0.844 

-1.147 
2.651 

-0.442 



TABLE II: Hopping amplitudes between different atom and orbital 
types (all numbers are in eV units). 

Table H] shows that the normalized bare energies of the Hg 
orbitals are on average larger than those of their Cd counter- 
parts. The HgTe layers are therefore more positively charged 
than the CdTe layers, resulting in a potential difference that 
lowers the orbital energies in the HgTe layers. This effect is 
taken into account by introducing a periodic potential U (r) 
which is added to all bare energies at position r. The potential 
depends on z only, and we write it in the form 



U{z) = —Uq cos 



27r 



2 



(12) 



where di 2 = aNi,2 are the thicknesses of the HgTe and CdTe 
layers, and d = di +^2 is the periodicity of the superlattice. In 
the symmetric case when Sq ~ 0, the potential function U{z) 
reaches its minimum in the middle of the HgTe layer and its 
maximum in the middle of the CdTe layer. However, we as- 
sume a certain asymmetry in U{z) which is characterized by 
the displacement So of these extrema. When the multilayer 
is grown under reasonable experimental conditions, such an 
asymmetry is inadvertently present due to the specific growth 
direction. For example, it is possible that the HgTe and CdTe 
materials are more likely to form an alloy at one of their in- 
terfaces. Since the resulting asymmetry is probably small, we 
take < \6o\ < a in the rest of this section. Also, by compar- 
ing the bare energies in Table|T]we estimate that the amplitude 
of the potential is Uq ~ 0.1 eV. 



B. Linear four-band approximation 

The Hamiltonian presented in Eqs. ^ and ( fTOl i can be 
solved by exploiting translational invariance and introducing 
the corresponding momentum k — {k^, ky, kz). On the other 
hand, the large periodicity in the z direction means that the 
Hamiltonian is represented by a large M x M matrix where 
AI ~ 40{Ni + N2)- It can therefore only be solved numeri- 
cally, and for relatively small layer thicknesses Ni^2- 

However, despite the complexity in this model, some of its 
properties can be deduced by referring to symmetry only. In 
the case of [/q 7^ and 60 7^ 0, the basic symmetries of the 
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system are time reversal symmetry (T) and the reflection sym- 
metries (R1.2) across the {y' , z} and the {x' , z} planes. The 
natural coordinates a;' = {x + y) / and y' = (.t — y)/V2are 
introduced to make these symmetries more explicit. Note that 
the reflection symmetries Ri 2 also lead to a two-fold rota- 
tional symmetry (S) around the z axis. In terms of symmetry, 
the tight-binding model studied in this section is equivalent to 
the superlattice model in Section HTPl 

The presence of both T and S symmetries puts two cru- 
cial restrictions on the band structure. First, the states on the 
kx' = ky> = line are two-fold degenerate. Since the num- 
ber of occupied bands is always even, the highest occupied 
band (HOB) and the lowest unoccupied band (LUB) have dif- 
ferent energies on this line with the Fermi level lying between 
them. Second, the band structure is invariant under the reflec- 
tion kz ^ — fcz. As shown in Section ITlI CI this implies that 
robust band touching between the HOB and the LUB is only 
possible in the k^ ~ and the kz = Tr/rf planes. 

Numerical investigation of the model indicates that band 
touching between the HOB and the LUB always occurs close 
to the kx' = ky/ = line. In perspective of this and the 
symmetry considerations above, we introduce simplified four- 
band models around the two special points at A; = (0,0, 0) and 
k = (0, 0, Tr/d). We only keep the nearest two bands on each 
side of the Fermi level, and assume that they are linear in the 
relative momentum Sk = {kx/, ky>,6kz) with respect to the 
corresponding special point. More formally, we project the 
Hamiltonian H onto a subspace spanned by four basis states: 
the appropriate eigenvectors of the full model at 5k = 0. The 
reduced Hamiltonian H is then a 4 x 4 diagonal matrix at 
6k = 0, and the linearity of the band dispersion is achieved 
by additional terms that are linear in 6k. 

At the special point, the reduced Hamiltonian can be written 
as Ho = -£(^^(1) + £;(2)r(2) where r^i'^) = (1 ± t^)/2 and 
£:(i,2) fjjg energies of the LUB and the HOB at 6k = 0. 
The energy levels are pairwise degenerate, and this degener- 
acy is split by a finite k^' or ky/ but not by a finite 6kz ■ To rep- 
resent this splitting, we need to add coupling terms between 
states corresponding to the same energy at 6k = 0. In the most 
general case, the additional terms in the Hamiltonian read 
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Ri 



R2 



1^1 
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(13) 



where the coefficients ax,y can be obtained from a compari- 
son with the full model. The choice of the Pauli matrices a^-^ 
corresponds to defining the basis states within the degenerate 
subspaces of Ha in a particular way. 

There is also coupling between states corresponding to dif- 
ferent energies at 6k = 0. Since the Hamiltonian H must be 
invariant under all symmetry operations of the system, there 
are only a small number of such coupling terms allowed by 
symmetry. The parities of the possible terms under the sym- 
metry operations are summarized in Table |llll Note that the 
parities of and cr^'^ are determined by the already estab- 
lished terms Ho and Hs which must be even under all sym- 
metry operations. Furthermore, the relations between differ- 
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kyl 



Sk, 
a 



TABLE III: Potential terms in the Hamiltonian and their parities un- 
der the symmetry operations of the system: time reversal (T), reflec- 
tion across the {y', z\ plane (Ri), reflection across the {x', 2} plane 
(R2), and two-fold rotation around the z axis (S). 



ent Pauli matrices imply that we only need to choose the par- 
ities of under T and Ri 2- The choice of these parities 
corresponds to setting the relative complex phases of the ba- 
sis states. Under the choice presented in Table HIH the most 
general contribution to the Hamiltonian takes the form 



Ud = {^xky'G^ + ^ykx'oy + ^z6kza^\ 



(14) 



where the coefficients fix,y,z are again to be determined from 
a comparison with the full model. The reduced Hamiltonian 
finally reads H = Hq + Hs + ^.0- It is a considerable simpli- 
fication with respect to H, and it only contains 9 parameters 
that need to be extracted from the full model. 



C. Conditions for robust band touching 

Band touching between the HOB and the LUB occurs in the 
full model when the two middle eigenvalues are equal in the 
simplified model. It can be shown that for a 4 x 4 matrix of the 
form H, this is possible if and only if the direction of the vec- 
tor B ~ {Pxky' T Pykx' , fiz6kz) lies halfway between the di- 
rections of the vectors A^-*^'^) = [a'x''^^ kx' joiy'^'^ kyl The 
two bands then cross each other as \6k\ is increased without 
changing the direction of 6k, whereas anti-crossing happens 
otherwise. Since the above condition requires the three vec- 
tors to lie in the same plane, the third component of B has to 
vanish. Due to /3z 7^ in general, we find that robust band 
touching can only occur in the 6kz ~ plane. 

Restricting our attention to this plane simplifies the problem 
because and B become 2D vectors. If we change the 

ratio kyl /kx' gradually from to 00, the ratios of the corre- 
sponding components in A^^-'^^ change in the same direction, 
while those in B change in the opposite direction between 
and ±00. This means that whether band touching happens at 
any 6k \s determined entirely by the signs of the different pa- 
rameters. Since we always choose ax,y > 0, the condition 
becomes straightforward: band touching occurs if and only if 
Px and Py have the same sign. 
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Let us now consider the special case of the symmetric po- 
tential with Sq = 0. By repeating the symmetry considera- 
tions in Section HUB I and taking into account the additional 
four-fold roto-reflection symmetry around the z axis, we find 
that the parameters from the full model are no longer indepen- 
dent because ai^'^-* = o^y''^^ and = ^Py This shows that 
band touching can only occur in this scenario if at least one of 
these parameters vanishes. However, the corresponding band 
touching is not robust because it requires the fine-tuning of 
a parameter. We conclude that robust band touching requires 
the asymmetry characterized by 5o 7^ 0, and expect that it 
becomes easier to observe as Uq and Sq increase. 



D. The Weyl semimetal phase 

The detailed behavior of the system is determined by how 
the coefficients from the full model depend on the external pa- 
rameters. Since this dependence is affected by the complexity 
of the full model, its understanding requires a numerical treat- 
ment. In perspective of this, we numerically investigate the 
phenomenon of robust band touching in the function of the 
layer thicknesses iVi 2, the strain eo = e(HgTe) in the multi- 
layer structure, the amplitude Uq of the superlattice potential, 
and the asymmetric displacement Sq. 

We first consider the dependence on the strain. If Uq ^ 
and So ^ 0, there are two ranges in eo close to zero with band 
touching in the A;^ = and the ^ Ti/d planes, respectively. 
The corresponding band touching is robust because it remains 
intact for an infinitesimal change in any of the external param- 
eters eo, Uq, and 5q. The upper and lower limits of the ranges 
are functions of C/q and 5^ as illustrated in Fig. |2] and we ver- 
ify the expectation from Section lTlICI that the ranges increase 
with both Uq and 5q. For the reasonable values of Uq ^ 0.1 
eV and 5q ^ a/2, the ranges are Aeo ~ 0.002. 
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FIG. 2: Critical strains eo against Uo at constant So = a/2 (left) and 
against So at constant Uo = 0.1 eV (right). The phase boundaries 
separate four distinct phases: the normal insulator (NI), the strong 
(3D) topological insulator (TI), the weak (2D) topological insulator 
(WI), and the Weyl semimetal (WS). The layer thicknesses are con- 
stant A^i = 3 and = 4 in both subfigures. 



Now we turn our attention to the layer thicknesses. Keep- 
ing the HgTe thickness A^i = 3 constant and varying the CdTe 
thickness N2 between 4 and 8 reveals that an increase in N2 
decreases Aeo. This is intuitive because 8q becomes smaller 
in comparison to d. Keeping the CdTe thickness N2 = ^ 
constant and varying the HgTe thickness A^i between 3 and 
7 shows that an increase in A^i shifts the ranges in eo down- 
wards. This means that the phases with robust band touching 
appear at more negative strains. 

To conclude that these phases are indeed Weyl semimetals, 
they need to satisfy one more condition: the lack of band over- 
lap. Even if there is robust band touching between the HOB 
and the LUB, the band structure becomes metallic if the high- 
est overall energy of the HOB is larger than the lowest overall 
energy of the LUB. It is an empirical observation that the indi- 
vidual band structures of the fcz = and the kz = ir/d planes 
are metallic when eo is sufficiently negative. This occurs for 
all eo < in the limit of C/o ^ or (5o — 0, while the crit- 
ical eo becomes slightly negative at larger values of Uq and 
5q . Furthermore, the appropriate bands of the k^ ~ and the 
kz = T^/d planes can overlap with each other as well. Since 
the band touching energies are different in the two planes, this 
typically occurs when there is band touching in one of the 
planes and almost band touching in the other one. 

We are now in the position to discuss the other phases 
around the Weyl semimetals. The overall transition between 
the two bulk phases at small and large values of eo is a 
2D topological phase transition because it happens via band 
touching around both special points on the axis. This 
means that the bulk phases on the two sides of this transi- 
tion do not have a 3D topological character: we identify them 
as the NI phase and the weak (2D) topological insulator (WI) 
phase. Since the spin-orbit coupling is stronger in HgTe than 
in CdTe, we argue that the system is in the NI phase when 
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FIG. 3: Phase diagram of the system against the strain eo and the 
HgTe thickness A'^i. The other parameters are constant: Uo = 0.2 
eV, Sq — a/2, and ~ 4. The phase boundaries separate five 
distinct phases: the normal insulator (NI), the strong (3D) topological 
insulator (TI), the weak (2D) topological insulator (WI), the band 
overlap metal (M), and the Weyl semimetal (WS). The dashed lines 
indicate approximate phase boundaries. 
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the HgTe layers are thin and in the WI phase when the HgTe 
layers are thickfi On the other hand, the phase between the 
Weyl semimetals is related to each bulk phase by a topologi- 
cal phase transition that has a 3D character because it happens 
via band touching around only one special point. We con- 
clude that this phase in the middle is the TI phase. Note that 
the NI and the WI phases are equivalent in terms of their 3D 
topology, therefore there is no need to distinguish between the 
two Weyl semimetals: they are intermediates in two equiva- 
lent phase transitions. 

The phase diagram of the system against the strain eq ™d 
the HgTe thickness iVi is presented in Fig. [3] Since its bound- 
aries are interpolated from only 5 points corresponding to in- 
teger values of Ni, the phase diagram is only correct on the 
qualitative level. Nevertheless, it provides useful guidelines 
for the realization of the Weyl semimetal phase in this multi- 
layer structure. The strain eo has to be positive to avoid band 
overlap but not too large because that would be hard to achieve 
experimentally. This gives a restriction on the thickness of the 
HgTe layers: the ideal dimensionless thickness of 4 < iVi < 6 
corresponds to an actual thickness of c?i ~ 3 nm, which is on 
the border of experimental reasonability. 



E. Connection witli the superlattice model 



IV. PHYSICAL CHARACTERISTICS 
A. Conductivity anisotropy 

It was shown in Ref. [H that the Weyl semimetal phase is 
metallic: when impurities are present, its conductivity is a fi- 
nite constant in the limit of zero temperature. Using the Boltz- 
mann equation, one finds a conductivity a = e'^v'^ /Qttj for 
each Weyl point, where v is the effective Fermi velocity and 
7 is the strength of the impurity potential. This finite con- 
ductivity is a characteristic experimental feature, especially in 
contrast with the neighboring NI and TI phases. In this sub- 
section, we demonstrate that the finite conductivity at T 
becomes highly anisotropic when the Weyl semimetal phase 
occurs due to broken inversion symmetry. 
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To illustrate the relationship with the results obtained in 
Section ini we discuss the arrangement of the Weyl points in 
the Weyl semimetal phase. There are 4 Weyl points that are 
related to each other by the symmetries of the system. As eo is 
gradually increased, and the transition from the NI (TI) phase 
to the TI (WI) phase happens through a Weyl semimetal, the 
Weyl points first appear at the k^' = line, move on approx- 
imately circular curves, and finally disappear at the ky' = 
line. This is in perfect agreement with the corresponding ar- 
rangement for the superlattice model in Section HTD] Indeed, 
the two models presented in Sections [II Dl and Fill Al obev the 
same symmetries, therefore it is understandable that their low- 
energy features are equivalent. 

The comparison of the band structures in Sections [II Al and 
IIIIBI also makes it possible to estimate reasonable values for 
the superlattice parameters in Eq. ([T]). The inversion-breaking 
voltage V corresponds to the energy difference - £'(2) 
between the HOB and the LUB, which is typically about 
0.05 eV in the tight-binding model. The Fermi velocity 

of the surface states becomes vp ~ c^i^if^ ~ 10^ ms^^, 
and the hopping amplitudes are estimated from the typical 
energy scale along the k^ axis: At^n ^ 10~^ eV. The 
small magnitude of At,n indicates that the band structure is 
relatively flat in the kz direction. Since Vr = vp ^ 10^ ms^^, 
and 10 nm gives = dy/ApAN ^ 10^ ms~^, this 
results in the relation Vz ^ v,- between the effective Fermi 
velocities around the Weyl points. Note that the energy scale 
At,n also translates into a maximal temperature T 10 K at 
which the Weyl semimetal phase is experimentally observable 
in this multilayer structure. 



FIG. 4: (Color online) Variation in the principal conductivities axx 
(red solid line) and ayy (blue dashed line) during a transition between 
the normal insulator (NI) and the topological insulator (TI) phases 
through the Weyl semimetal phase (WS). The transition parameter is 



i;|.(A*^' - AI^')/[V^{St - Sn)] and the conductivities ai'e 



measured in units of ctq = 2e vp/STvy. 



To achieve this, we consider the model in Section FTTPl and 
derive an expression for the conductivity tensor in the limit of 
small St — S^. When the condition of the Weyl semimetal 
phase is satisfied, there are 4 Weyl points at angles 9i — 6, 
02 = —0, 03 = TT + 0, and 04 = TT — 0. Due to the convention 
< < tt/2 we find that gradually decreases from 7r/2 to 
during a transition from the NI phase to the TI phase. For 
each Weyl point labeled by /, the conductivity tensor in the 
{x, y, z) basis takes the form^ 



<i 



v"^ cos 01 sin 01 



67r7 



cos Ui sm Ui 




vl SIT? 01 









(15) 



where we exploit vt ^ Wr relating the effective Fermi veloci- 
ties. Adding the contributions of all 4 Weyl points, there is a 
cancelation in the off-diagonal terms, and we obtain 



4 ^ ■? vi COS. 

a, = — vlsm 

'=1 V 



(16) 



As the transition between the NI and the TI phases takes place 
through the Weyl semimetal phase, the conductivities in the 
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X and y directions change in a complementary fashion. In 
particular, aj.j. vanishes on the Nl side and ayy vanishes on the 
Tl side of the Weyl semimetal phase. For an illustration, see 
Fig. |4] The conductivity in the z direction is approximately 
constant with azz <C crxxT<^yy due to ^ Vr- Such a strong 
conductivity anisotropy that depends sensitively on the system 
parameters is a potential hallmark of a Weyl semimetal with 
broken inversion symmetry. 



sponding dimensionless energies E = ±-\/k^ + are in- 
dicative of surface states with Dirac dispersion between NI 
and TI phases of matter. In the more relevant case ofV^Q, 
these analytic solutions only find straightforward generaliza- 
tions for K,. = when 



I*) = (i,T«,Tl,l)exp( -^ii^y 



(20) 



B. Topological surface states 

Since Weyl semimetals are topological phases of matter, 
they are characterized by topological surface statesi^ In this 
subsection, we consider the model in Section HTD] and demon- 
strate the existence of these surface states. Although we 
choose a specific situation and also make a couple of simpli- 
fying assumptions in the following, the topological nature of 
the surface states ensures that they exist under more generic 
circumstances as well. 

In our specific situation, the interface is in the {x, z} plane, 
therefore any spatial variation is in the y direction only. This 
implies that and are still valid quantum numbers. Since 
Ay TV in Eq. (|7]) do not depend on ky, we can determine 
the surface states without taking the explicit k dependence 
into account, and then simply substitute the appropriate val- 
ues /S.T,N for each k^- It is assumed that only At changes 
with y and the other parameters are constant: Ay < Ajv at 
y -> — oo. At = Ajv at ?/ = 0, and Ay > Aat at y — >• +oo. 
Furthermore, if At changes sufficiently slowly, we can ap- 
proximate it with a linear function in the important region 
around i/ = 0: we write At — A^r = Ky. Expanding the 
kz dependent terms up to first order in k'^ = /c^ — tt/ d, we 
find that the surface states \^) with energy E need to satisfy 



E\^) 



VFT'{-icj''dy~aykx) + VT'' 



KyT"" + k'^dA 



NT" 



(17) 



along with l^*) in the limits of y ^ ±oo. To make the 
subsequent discussion of the surface states more transparent, 
we introduce the dimensionless form 



(18) 



where the variables y 

Kx ~ vpkx/A, and k 
while A 



y^K/vF. E = E/k, V = 1//A, 
= k'^dApi/A are all dimensionless, 
\/Kvf is a characteristic energy scale. 
As a starting point in our discussion, we consider the limit 
of V = 0. In this case, there are two distinct surface state 
solutions for each Kx and that take the analytic form 



I*) = {^,^ie'^,^e'^,l) exp 



y 



(19) 



in the (T^, TJ,, B^, BJ,) basis, where ta.nLp ~ Kx/kz and 
the letters T/B stand for the top/bottom surfaces. The corre- 



and the dimensionless energies are E = ±Kz- Note that these 
surface states decay in an oscillating fashion at y — > ±oo, and 
the wave vector ky = ±V/vf of the oscillations corresponds 
to the radius of the circle in the kz — Tr/d plane along which 
band touching occurs in SectionlU 

In the most generic case of V ^ and Kx ^ 0, we solve 
Eq. ( fTSl l numerically and find that there are still two distinct 
surface states j'l') for each Kx and Kz- The ratios of the vec- 
tor components in \'^) are no longer independent of y, which 
explains why simple analytic solutions like those in Eqs. ( fT9] l 
and ( |20] | can not be obtained. We verify that the surface states 
follow a Dirac dispersion at small momenta Kx^z ^ 1, even 
when the dimensionless voltage V is large. However, the ef- 
fective Fermi velocity in the Kx direction is reduced by a fac- 
tor that is empirically found to be exp{—V^), and hence the 
dispersion relation at small Kx^z becomes 



(21) 



Unlike in the direction where the analytic solution guaran- 
tees the linearity of the dispersion for all Kz, there is a devi- 
ation from the linear dispersion in the Kx direction. As illus- 
trated in Fig. |5] the Dirac dispersion in Eq. (|2T]) is only valid 
for small enough Kx- 
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FIG. 5: Dimensionless energies of the surface states in the function 
of the momentum k,x when k,z = Q and V = 1. The dashed lines are 
linear asymptotes given by Eq. i2H in the Kx <^ 1 limit. 

If we now restore the k dependence of At.w, the surface 
states 1 5') remain the same for each kx- However, they only 
exist at those k^ for which the difference At — Ajv changes 
sign between y ±oo. Assuming without loss of generality 



that St > Sn in Eq. ^ and that A^' — A^'' does change 
sign, we find that At > A^v for all kx al y ^ +oo, while 
At < Atv is only true aty ^ —oo for |/ca; | < fco- The critical 
momentum ko marks the equality At = A^y atj/ — oo. 
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which is one of the conditions required for band touching in 
Section While the material at y > is definitely in the 
NI phase, the material at y < is in the Weyl semimetal 
phase if fco < V/vp so that band touching occurs, and it is 
in the TI phase if fco > V/vp so that band touching does not 
occur In the former case, surface states exist between the 
coordinates = ±fco of the Weyl points. As the Weyl points 
first appear at the ^ line, and then start to move further 
away from it, the range in k^ increases and more surface states 
appear. Remarkably, this range characterized by fco further 
grows when the material at y < is akeady in the TI phase, 
and there is no band touching at all. 

V. SUMMARY 

We proposed a time reversal invariant reaUzation of the 
Weyl semimetal phase that occurs due to broken inversion 
symmetry. We considered both a superlattice model adapted 
from Ref. and a tight-binding model describing an experi- 
mentally reasonable HgTe/CdTe multilayer structure. The su- 
perlattice model was suitable for analytic calculations due to 
its simplicity, while the more realistic tight-binding model re- 
quired a numerical treatment. 

Although the formulations of the two models are very 
different, their identical symmetries lead to equivalent low- 
energy features. It should be remarked that, as seen from 
the generality of the superlattice model, the Weyl semimetal 
could be achieved in many possible material structures. 
Exploration of potential compounds other than HgTe/CdTe 
would be extremely interesting, especially given the need to 
tune strain in the latter to observe the desired physics. 



For both models considered, we found a Weyl semimetal 
phase between the NI and the TI phases. This phase is char- 
acterized by a finite number of Weyl points with robust band 
touching at the Fermi level: the band touching occurs for a 
finite range of the system parameters, and hence it can not be 
removed by an infinitesimal perturbation. We further verified 
that the band touching points are proper Weyl points with a 
linear dispersion relation around them. 

In terms of experimental observation, the potential hall- 
marks of the Weyl semimetal phase with broken inversion 
symmetry are a strong conductivity anisotropy and the 
presence of topological surface states. The highly uncon- 
ventional low-temperature and low-frequency bulk transport 
is discussed in Ref. |3 The Dirac dispersion relation of 
the surface states is indicative of TI materials, but these 
surface states only exist in a region of momentum space that 
is determined by the positions of the Weyl points. The Weyl 
semimetal phase between the NI and the TI phases described 
in this paper is therefore qualitatively new in terms of its 
topological surface states as well. 

Acknowledgments 

We are grateful to T. L. Hughes for useful discussions. This 
research was supported in part by NSF grants DMR-0804564 
and PHY05-5 1 164 , and by the Army Research Office through 
MURI grant No. W91 l-NF-09- 1-0398. G. B. H. acknowl- 
edges the support of J. DriscoU (Trinity College) and the hos- 
pitality of KITP during this work. 



' M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). 
^ D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, 

Phys. Rev. Lett. 49, 405 (1982). 
^ C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005); B. 

A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 106802 (2006). 
" B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 314, 1757 

(2006) . 

^ J. E. Moore and L. Balents, Phys. Rev. B 75, 121306(R) (2007); 
L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 106803 

(2007) ; R. Roy, Phys. Rev. B 79, 195322 (2009). 

^ M. Konig, S. Wiedmann, C. Brune, A. Roth, H. Buhmann, L. W. 
Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 318, 766 (2007); 
D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, and M. 
Z. Hasan, Nature (London) 452, 970 (2008). 



^ X. Wan, A. M. Turner, A. Vishwanath, and S. Y. Savrasov, Phys. 
Rev. B 83, 205101 (2011); K.-Y. Yang, Y.-M. Lu, and Y. Ran, 
Phys. Rev. B 84, 075129 (2011). 

A. A. Burkov and L. Balents, Phys. Rev. Lett. 107, 127205 (201 1). 
' S. Murakami, New. J. Phys. 9, 356 (2007). 

C. Herring, Phys. Rev. 52, 365 (1937). 
" X. Dai, T. L. Hughes, X.-L. Qi, Z. Fang, and S.-C. Zhang, Phys. 

Rev. B 77, 125319 (2008). 

A. Kobayashi, O. F. Sankey, and J. D. Dow, Phys. Rev. B 25, 6367 
(1982). 

J. N. Schulman and Y.-C. Chang, Phys. Rev. B 33, 2594 (1986). 
A. A. Burkov, M. D. Hook, and L. Balents, Phys. Rev. B 84, 
235126 (2011). 



